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1. Introduction 


Gait analysis of patients with neurological disorders, including multiple sclerosis (MS), is impor¬ 
tant for rehabilitation and treatment. In a clinical setting, gait analysis performed by physicians or 
therapists involves observing a patient’s gait and manual documentation of subjective assessments. 
Different clinical scores are proposed in order to quantify the level of progression of MS. Among 
them the expanded disability status scale (EDSS) and the self-scored multiple sclerosis walking scale 
(MSWS) are widely used in clinical practice pLl [ 2 j . However, these clinical scores are subjective and 
are unable to detect subtle changes in a subject’s gait due to the progression of the disease or the 
response to treatment. A more systematic gait analysis can be carried out in a gait laboratory 
using motion capture systems, force-plates, and electromyography (EMG) sensors. However, the 
required set up involves elaborate preparation and marker placement, rendering it unsuitable for 
use as a point-of-care technology. 

A number of indices based on time or distance characteristics of human gait cycle have been 
proposed for objective gait assessment in neurological patients with gait abnormality. Specifically, 
studies have shown that a shorter stride length and gait swing time and a higher double support 
percentage in a gait cycle (i.e., the fraction of the time in a gait cycle where the two feet are on the 
ground) are observed in neurological patients with gait abnormality [SIS]- In addition, the range 
of the hip and knee angle is smaller in MS patients as compared to healthy individuals as reported 
in [7[l8]. Correlations between time-distance and joint angle indices and the EDSS score have been 
investigated as well um- 

Angles of the lower extremity joints in a MS patient’s gait have also been investigated. However, 
researchers have reported contradictory results. Specifically, reduced hip and knee flexion and ankle 
plantarflexion at heel strike (i.e., the point in the gait cycle when foot reaches the ground), and hip 
and knee extension at toe-off (i.e., the point in the gait cycle where foot is no longer in contact with 
the ground) are reported in [SJE!], while increased hip and knee flexion and ankle plantarflexion at 
heel strike has been observed in [8]. The analysis in these studies involve point-to-point comparisons 
between preselected peak points of the joint angles during a gait cycle. 

Although accurate gait analysis for MS patients could provide valuable insight on a patient’s 
condition and the progression of the disease, an easy-to-use system which can encourage clinical 
adoption does not exist. Motion capture systems, which are mainly used for experimental data 
collection for gait analysis, are expensive. Furthermore, motion capture systems can only track a 
set of reflective markers mounted on certain anatomical landmarks of a patient’s body, and hence, its 
use requires a time-consuming patient preparation process. Moreover, operator training for motion 
capture systems is necessary to ensure acceptable accuracy. Finally, these systems are not portable 
and require a dedicated space which makes them further impractical for clinical applications. 

The Mircrosoft Kinect® sensor, which was developed for motion recognition in gaming appli¬ 
cations, is an ideal candidate for an inexpensive system providing the capability for human gait 
analysis. The Kinect sensor includes a color camera and a depth sensor, consisting of an infrared 
projector and camera, and provides full-body 3D motion capture. The Kinect sensor has been 
used for various clinical and non-clinical applications. The authors in m use the Kinect sensor for 
pose identification. In |12] , joint angles identified by the Kinect sensor are compared to the “gold 
standard” obtained by a marker-based motion capture systems for healthy subjects, which showed 
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reasonable accuracy for clinical applications. The validity of the Kinect sensor for the assessment 
of the postural control was examined by comparing the result with the result of a marker-based 
motion capture system too m- Also, the accuracy of this sensor for movement measurement in 
people with neurological disease, such as Parkinson, has been examined nu 

However, a limited number of studies have been performed to investigate the feasibility of Kinect 
specifically for gait analysis of MS patients |15Ul7j . In [17] . the short maximum speed walk test 
was proposed by the authors to be measured with Kinect, where the correlation of this index with 
EDSS was investigated. Furthermore, based on machine learning and image processing techniques, 
movements of MS patients were compared with healthy subjects to identify subgroups with similar 
movement patterns [16]. Authors in [18] have developed a framework to identify MS patients from 
healthy subjects by analyzing a number of tasks such as finger-to-nose and finger-to-finger tests. 

In this research, we develop a framework to quantify the gait abnormality of MS patients using 
a Kinect for Windows (version 1) camera. We show that previously introduced indices quantifying 
gait abnormality as well as an index for quantifying a patient’s gait pattern introduced in this report 
are in strong agreement with clinical scoring systems that quantify the degree of the progression of 
the disease. Specifically, we show that there is a correlation between a number of gait indices and 
the self-scored MSWS as well as the clinical ambulation score mm when clinical data is obtained 
using a Kinect camera. The clinical ambulation score, which is one of the required scores for 
computing the EDSS score, is determined by a clinician through clinical observation of a patient’s 
gait. Considering the clinical ambulation score appears to be more meaningful for patient gait 
abnormality assessment compared to EDSS which quantifies a patient’s general disability due to 
the MS disease and involves factors other than gait. 

In this report, in addition to the previously introduced gait indices, a novel set of MS gait 
indices based on the concept of dynamic time warping |T9] is introduced. The newly introduced 
indices can characterize a patient’s gait pattern as a whole (rather than considering isolated events 
in a gait cycle), and quantify a subject’s gait “distance” from the healthy population. We will 
investigate the correlation of these novel indices with the MSWS and the clinical ambulation score. 
This work establishes the feasibility of using the Kinect sensor for clinical gait assessment for MS 
patients. 


2. Gait Analysis for MS Patients 

2.1. Subjects 

In this study, 10 male and female MS patients (9 females and 1 male), and 10 sex and age 
matched healthy control subjects were asked to walk in front of the Kinect camera for 5-10 trials 
(i.e., a video sequence which involves a subject moving in front of the camera on a straight line). The 
best five captures were selected for the analysis. This study was approved by the Research Ethics 
Board at the McGill University Health Center and was conducted at the Montreal Neurological 
Institute. For each patient the MSWS and the clinical ambulation score were assigned by an on-site 
physician on the day of the study. MS patient information and their clinical assessment scores are 
summarized in Table [TJ whereas Table [2] summarizes the information related to healthy control 
subjects. 
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Table 1: MS Subjects Specifications 


Patient No. Sex 

Age 

Ambulation Score 

MSWS 

Height [cm] 

Weight [kg] 

PI 

F 

53 

6 


89.6 

146 

55 

P2 

F 

41 

0 


25 

170 

86 

P3 

M 

79 

5 


47.9 

169 

91 

P4 

F 

69 

1 


60.4 

150 

81 

P5 

F 

75 

5 


72.9 

157 

89 

P6 

F 

60 

6 


81 

170 

79 

P7 

F 

55 

4 


83.3 

168 

81 

P8 

F 

70 

9 


97.9 

178 

81 

P9 

F 

53 

6 


75 

174 

59 

P10 

F 

55 

1 


45.8 

161 

51 



Table 2: 

Healthy Control Subjects Specifications 



Control Subject No. Sex 

Age 

Height [cm 

] Weight [kg] 




Cl 

F 

53 

160 

81 




C2 

F 

36 

165 

70 




C3 

M 

80 

185 

135 




C4 

F 

71 

163 

82 




C5 

F 

72 

159 

50 




C6 

F 

62 

147 

59 




C7 

F 

57 

170 

63 




C8 

F 

67 

170 

63 




C9 

F 

50 

168 

69 




CIO 

F 

51 

165 

53 



In each trial, subjects were asked to walk at a normal pace. They were instructed to start their 
gait outside of the camera’s field of view in order to ensure that their normal walking patterns have 
been established once they reach the capture zone. As the Kinect field of view is limited, depending 
on the stride length of the subjects, one or multiple gait cycles might have been captured. The 
stored information related to each subject was de-identified to ensure patient privacy. Depending 
on the level of progression of the disease, patients may have walked with or without the use of 
assistive devices. Normal subjects passed an interview with the on-site physician to ensure they 
did not suffer from any gait abnormalities and do not possess any condition that could affect their 
gait. 
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2.2. Data Analysis 

The Kinect for Windows sensor with the use of its software development kit (SDK) provides 
three-dimensional skeletal data on 20 joint positions over time. For the lower extremity, these joints 
consist of pelvis, knee, and ankle for each leg. Kinect captures the video up to 30 frames per second. 
Joint positions are expressed in an inertial reference frame in which the y-axis is in the direction 
of the runway, the z-axis is perpendicular to the ground, and the x-axis is mutually perpendicular 
to both (see Figured]). 



Figure 1: Captured data and the identified joints using Microsoft Kinect. The left hip angle 0 lh and left 
knee angle 0 lk are shown in the Figure as well. 


The subject’s kinematic properties can be extracted using the value of the joint positions. Joint 
angles, and if required, angular velocities and accelerations can be calculated based on the time- 
history of the joint positions. The first step in the process involves identifying a complete gait cycle 
in the captured data. This can be accomplished by considering ankle position variations over time. 
A gait cycle starts with heel strike, implying that the ankle joint position is stationary (see Ths in 
Figured])- At toe-off the leg starts to swing, and hence, the ankle position starts to change (see Tto 
in Figured])- Finally, the gait cycle terminates by the terminal swing in which the ankle position 
comes to rest again (see Tts in Figured]). The ankle joint position variation for a representative 
subject is shown in Figured] 

In order to distinguish an abnormal gait from normal, and potentially quantify its degree of 
abnormality, appropriate gait indices need to be defined. Three categories of indices are defined for 
gait analysis in this study; namely, time-distance indices, joint angle, which have been discussed 
in the literature, and gait pattern, which is introduced in this report. Here, four time-distance 
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Figure 2: Ankle joint position variation captured 
for Patient 6 (Trial 5). Heel strike, toe-off and 
terminal swing phases are denoted by Ths, Tto, 
and Tts, respectively. 



Figure 3: Hip angle variations for an MS patient 
and a healthy control subject. 


indices, namely, subject velocity, stride length, stance percentage, and stride width are considered 
as suggested by authors in [Mi- In the joint angle category, hip and knee range of motion is 
considered due to the fact that studies show that MS can substantially affect these joints [7]. 
Finally, a novel index, which can capture the general gait pattern of the subject and its deviation 
from a healthy gait, is introduced. 

A subject’s velocity for one gait cycle can be calculated by dividing the stride length to the 
gait cycle time (from Ths to Tts)- To be able to meaningfully compare distance indices, they 
are normalized based on a subject’s height to compute normalized velocity V n . Furthermore, 
normalized stride length L n , defined as the distance travelled by the ankle in one gait cycle (i.e., 
between the heel strike Ths and the terminal swing phases Tts), stride width W, defined as the 
distance between the two ankles in the double support phase projected on the x-axis, and stance 
percentage S, defined as the ratio of the stance time to the gait cycle time, can be calculated. 
Finally, the knee and hip angles need to be computed using the location of the joints. Knee angle 
is defined as the angle between the thigh and the leg segments, and the hip angle is defined as the 
angle between the z-axis and the thigh segment (see Figured]). 


2.2.1. Dynamic Time Warping 

Although minimum/maximum joint angle values, range of motion, and time-distance indices 
can provide valuable information on a subject’s gait characteristics, they only provide a snapshot 
at a specific time instant. Developing a framework capable of analyzing a complete gait cycle (as 
opposed to analyzing certain points in a gait cycle) provides a holistic approach for gait analysis 
and can complement the information provided by other indices. This involves developing a distance 
metric to quantify the level of abnormality of a gait cycle with respect to healthy individuals. There 
have been prior attempts in comparing MS gait patterns with control subjects by point-to-point 
comparisons at certain gait phases such as heel strike or toe-off ISflD]. However, no study has 
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introduced a mechanism to compare complete gait cycles. 

Here, we propose a set of novel MS gait indices based on the dynamic time warping (DTW) 
framework. DTW, which was initially proposed for speech recognition applications [[20], provides a 
framework to find an optimal alignment between two time series that have different time scales. This 
is ideal for comparing sequences representing the human gait cycle as any gait cycle includes the 
same gait phases (i.e., heel strike, toe-off, etc.), however, the transition time between these phases 
varies from subject to subject. DTW has been previously used in other non-clinical contexts 
such as human motion recognition and in identifying different modes of movements and motion 
patterns [21]l22j. DTW defines a cost function and uses a nonlinear transformation to warp the 
two sequences in order to minimize the cost function. The optimal value of the cost function can 
be regarded as a “distance measure” between the two sequences. 

Let us consider two sequences A = {a\, 02 ,..., a]y}, o n £l,n = l, ..., N, and B = { 61 , bo,..., 
&m}> b m G M, m = 1, ..., M. A local distance cost between two elements of the sequence a n £ A 
and b m € B is a mapping c : A x B —>• M + such that c[a n ,b m ) increases as the mismatch between 
a n and b m increases. Next, we define a cost matrix C = [c nm \ € H NxM , where c nm = c(a n ,b m ), 
n = 1,... ,N, m = 1,... ,M. A sequence p = {p\, p 2 , ■ ■■, pl}, Pi = ( ni,mi ), l = 1, ..., L, 
ni € {1, ..., N}, mi € {1, ...,M} is referred to as a warping path if it satisfies the boundary, 
monotonicity, and step size conditions discussed in m- A warping path defines the point-to-point 
correspondence between the two sequences. The total cost c p (A,B) for the two sequences A and 
B along a warping path p is defined as 


L 

c p(A, B) .— ^ ^ c(o. n; , b mi ), a ni , b mi £ p. (1) 

1=1 

Consider the set of all possible warping paths for the two sequences A and B denoted by P. 
The dynamic time warping distance between A and B is defined 

DTW(A, B) := c p *(A, B), (2) 

where p* := argmin pg pC p (A, B) is the optimal warping path. 

2.2.2. Dynamic Time Warping for Gait Analysis 

Here, we investigate the application of DTW to quantify MS disease progression. Specifically, we 
investigate its correlations with the MSWS and clinical ambulation score. The DTW framework can 
be used to align two sequences, for instance unaligned hip angle variations over time for a patient 
and a normal subject gait. Figure [3] shows the the two sequences for a representative patient and 
a control subject. The value of the DTW distance defined in ([2]) shows the deviation of a patient’s 
joint pattern from a normal subject. A larger DTW indicates a larger deviation between patient 
and the healthy subject gait patterns, and potentially a more advanced stage of MS. In longitudinal 
studies, all of these novel DTW distances could potentially show the progression of the MS disease 
in a particular patient over time and provide insight on a patient’s response to treatments. 

In this research, in order to evaluate a patient’s gait, we extract time series for hip and knee 
angles for a complete gait cycle and compute the hip and knee DTW distances with respect to 
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a set of control subjects. Note that maximum hip and knee flexion angles are reported to be 
different in MS patients compared to healthy control subjects [BHlOj. In this study, the goal is 
to further investigate the relationship between the degree of gait abnormality and lower extremity 
joint variations over a complete gait cycle and not limit the analysis to extreme joint angles. 

The dataset included in this study is composed of gait data collected from n p patients, where 
each patient completed m p trials. In addition, the control dataset is composed of gait data col¬ 
lected from n c control subjects, where each control subject completed m c trials. As discussed in 
Section ^. 21 for each trial, one complete gait cycle from heel strike to toe-off is identified and the rest 
is discarded. Therefore, overall the dataset includes n p m p gait cycles corresponding MS patients 
and n c m c gait cycles corresponding to healthy control subjects. In our study, n p = n c = 10 and 
m p = m c = 5. 

For each patient, hip joint angles for the left and right legs over the entire gait cycle are calcu¬ 
lated and stored in arrays OhRij, 0 r,H; j, respectively, where i G {1 ,... , n p }, and j E {1 ,... , m p }. 
Similarly, knee joint angles for the patient’s left and right legs are stored as 0 lk ; ; #RKi respec¬ 
tively, where i € {1,... ,n p } and j € {1,... , m p }. For control subjects, hip and knee joint angles 
for the left and right legs are calculated and stored in t > &lk q r > ^rk,, r , respectively, 

where q € {1 ,..., n c } and r € {1,..., m c }. 

Here, we introduce an index referred to as the mean dynamic time warping distance. This 
index quantifies the degree of dissimilarity between a patient’s joint angle pattern compared to a 
set of control subjects. Specifically, the DTW distance between two time series, namely, a patient’s 
joint angle for a given trial and a set of control subjects joint angles for all available trials are 
computed. Note that these distances are computed for the left and right legs independently. That 
is, a patient’s left (right) leg joint angle sequence is compared with all left (right) joint angle 
sequences for all control subjects. Hence, overall n c m c DTW distances are computed for each leg 
(in our case n c m c = 50). 

Considering the fact that m p trials are available for each patient, m p n c m c = 250 DTW distances 
are obtained for each patient’s leg and each joint. Hence, a total of 2m p n c m p = 500 DTW distances 
are computed for each joint (i.e, knee and hip). These distances are then averaged to compute 
the mean dynamic time warping distance associated with the knee or hip joint for each patient. 
Specifically, for Patient i, i € {1 ,... ,n p }, the mean dynamic time warping distance for knee and 
hip joints, denoted by Dk p and Dh p , respectively, are defined as 


D 


K P 


Dh p := 


1 


m p n c m c 


EEE™ ( 8lK i J > 0LK,, ) + 


1 


m p n c m c 


m p n c m c ^ , 

1 = 1 (7=1 V 1 


EEE DTW ( 8 kk, j ,0 R k, 


m p n c m c “ 

1 1 = 1 <7=1 r =1 


1 


Trip n c m c 


EEEdt w (» IHu ,^„) + 


1 


n c m c 


m p n c m c 

L ^ j = l q= 1 r =1 


EEE d ™(«»» u *», 


m p n c m c 

^ j=l q =1 r =1 


( 3 ) 


( 4 ) 


Similarly, the mean dynamic time warping distance associated with the knee and hip joints for 
control subjects denoted by Dk c and Dh c , respectively, can be defined, where the distance between 
joint sequences for a control subject is compared with all other control subjects. 
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3. Results and Discussion 


The aforementioned time-distance and joint angle indices discussed in previous literature, and 
the newly introduced mean DTW indices associated with quantifying gait abnormality are com¬ 
puted for all patients and control subjects. A list of all the computed indices is given in Table [3j The 
mean value for each index (averaged over 5 trials) for patients and control subjects are summarized 
in Tables U and 0 respectively. 

The results are also illustrated in Figures HHTTl in form of box plots. The hip and knee posi¬ 
tions for Patient 9 showed severe artifacts, and hence, this patient was excluded from the analysis 
for indices involving hip and knee. Note that for MS patients, the median velocity and median 
stride length are smaller compared to healthy control subjects (Figures 0] and [5]), while the median 
stance percentage, as well as the knee and hip ranges of motions are larger (Figures [6][9]). These 
observations are in agreement with previously reported gait characteristics of MS patients M- 
Furthermore, the newly introduced indices based on the DTW distance can provide further insight 
on the general gait pattern. More specifically, the median hip and knee mean DTW distances in¬ 
crease significantly in subjects with MS disease (Figures 1 101 and Figure [TT]). Larger DTW distances 
imply that the patient’s gait pattern in the MS population is not “similar” to the healthy control 
group. 


Table 3: Gait Indices for MS Patients 


Index 

Symbol 

Unit 

Description 

Normalized velocity 

v n 

sec -1 

Velocity normalized by hight 

Normalized stride length 

L n 

Unitless 

Stride length normalized by hight 

Stance percentage 

S 

Unitless 

Stance time divided by gait cycle time 

Step width 

w 

m 

Distance between the two ankles in double support phase projected on the x-axis 

Hip range of motion 

«H 

deg 

Difference between minimum and maximum of the hip angle 

Knee range of motion 

a K 

deg 

Difference between minimum and maximum of the knee angle 

Mean dynamic time warping distance for knee 

D k 

deg 

See © 

Mean dynamic time warping distance for hip 

D h 

deg 

See © 


Table 4: Gait Indices for MS Patients 


Patient 

< 

3 

(X 

1 

i— 1 

Ln 

S[%] 

W[m] 

a: A'[deg] 

a A [deg] 

D k [deg] 

D h [deg] 

PI 

0.48 

0.44 

66 

0.68 

36.0 

17.9 

191 

182 

P2 

1.23 

0.78 

59 

0.61 

42.7 

25.8 

262 

156 

P3 

0.82 

0.67 

58 

0.59 

35.5 

22.0 

275 

178 

P4 

0.48 

0.51 

63 

0.81 

34.9 

24.1 

275 

207 

P5 

0.64 

0.58 

61 

0.84 

30.9 

23.1 

285 

155 

P6 

0.48 

0.55 

69 

1.04 

26.1 

20.0 

362 

272 

P7 

0.73 

0.67 

58 

1.02 

34.2 

26.7 

265 

178 

P8 

0.39 

0.47 

52 

1.01 

39.2 

19.5 

516 

369 

P9 

0.56 

0.59 

54 

0.83 

27.3 

32.1 

293 

178 

P10 

0.87 

0.63 

59 

0.42 

43.7 

25.7 

150 

142 


Next, we investigate the reliability of each index using the intraclass correlation coefficient 
(ICC). Specifically, the ICC is calculated for both the patient and control groups. The results are 








Table 5: Gait Indices for Control Subjects 


Control 

Vnfs- 1 ] 

L n 

S[%] 

W[m] 

a if [deg] 

ol h [deg] 

D k [deg] 

D H [deg] 

Cl 

1.35 

0.78 

45 

0.64 

42.9 

19.3 

161 

155 

C2 

1.12 

0.78 

52 

0.54 

50.1 

32.8 

168 

124 

C3 

1.12 

0.73 

53 

0.95 

33.8 

26.9 

206 

123 

C4 

0.97 

0.66 

53 

0.79 

37.8 

26.2 

323 

223 

C5 

1.38 

0.90 

48 

0.41 

50.5 

32.8 

160 

114 

C6 

1.25 

0.81 

49 

0.65 

42.6 

24.2 

145 

117 

C7 

1.19 

0.80 

48 

0.58 

52.7 

30.6 

164 

123 

C8 

1.21 

0.84 

46 

0.50 

46.6 

31.1 

156 

114 

C9 

1.16 

0.71 

54 

0.31 

43.5 

27.2 

195 

140 

CIO 

1.42 

0.77 

54 

0.62 

50.5 

38.1 

170 

147 


1.4 



o 

_o 

% 0.8 
N 

"co 

E 0.6- i 

o i 


0.4 


0.2 


Patients Controls 


Figure 4: Normalized velocity V n for pa¬ 
tients and control subjects. 
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Figure 5: Normalized stride length L n 
for patients and control subjects. 


given in Tables [6] and [7| As it can be seen in both groups all indices have acceptable intraclass 
correlations with the exception of step width in MS patients. Also, the knee range of motion index 
possess a wide confidence interval. In order to statistically compare the differences between the 
value of indices corresponding to the patient and control groups, the paired t -test was used (see 
Tables [6] and [71). All indices show statistically significant difference between the patient and control 
groups as signified by small p -values (p < 0.05). 

Finally, we note that all considered indices with the exception of stance percentage, and hip 
range of motion show a high degree of correlation with the clinical ambulation score and MSWS. 
The stride width shows an acceptable degree of correlation with MSWS but not with the ambulation 
score. This observations signify the fact that the Kinect sensor can provide an objective framework 
for assessing gait abnormality that can be of clinical utility. The data can provide invaluable insight 
on the progression of the disease and response to treatment. 
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Figure 6: Stance percentage S for pa¬ 
tients and control subjects. 
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Figure 8: Knee range of motion ax) for 
patients and control subjects. 
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Figure 7: Step width length W for the 
patients and control subjects. 
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Figure 9: Hip range of motion ckh) for 
patients and control subjects. 


Table 6: Statistics for time-distance-derived indices 


V„[s- 1 l 


MS Mean (SD) 

Control Mean (SD) 

MS ICC (95% Cl) 

Control ICC (95% Cl) 
Correlation Ambulation Score (95% Cl) 
Correlation with MSWS (95% Cl) 
t-test p -value 


0.4 (0.14) 

1.2 (0.14) 

0.99 (0.984,0.998) 
0.93 (0.84,0.98) 
-0.69 (-0.81,-0.50) 
-0.86 (-0.91,-0.76) 
io- 6 


L„ 

0.6 ( 0 . 1 ) 

0.8 (0.07) 
0.93 (0.83,98) 
0.90 (0.76,0.97) 
0.54 (0.30,0.71) 
0.69 (0.50,81) 
10~ 4 


S[%] 

0.6 (0.05) 

0.5 (0.03) 

0.71 (0.25,0.92) 
0.75 (0.38,0.93) 
-0.14 (-0.40,0.15) 
0.04 (-0.24,0.33) 
10~ 4 


W[m] 

0.8 ( 0 . 2 ) 

0.6 (0.18) 
0.50 (0.63,0.96) 
0.80 (0.52,0.95) 
0.43 (0.17,0.64) 
0.51 (0.26,0.70) 
0.044 
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Figure 10: Knee mean DTW distance 
Dk for patients and control subjects. 


Figure 11: Hip mean DTW distance Dn 
for patients and control subjects. 


Table 7: Statistics for angle-derived indices 



aft-[deg] 

<*H [deg] 

Dk [deg] 

D h [deg] 

MS Mean (SD) 

36 (5) 

23 (3) 

236 (68) 

210 (47) 

Control Mean (SD) 

47 (9) 

30 (4) 

191 (54) 

156 (36) 

MS ICC (95% Cl) 

0.61 (-0.01,0.90) 

0.92 (0.80,0.98) 

0.88 (0.69,0.97) 

0.78 (0.43,0.94) 

Control ICC (95% Cl) 

0.89 (0.74,0.97) 

0.98 (0.94,0.99) 

0.93 (0.82,0.97) 

0.82 (0.54,0.94) 

Correlation Ambulation Score (95% Cl) 

-0.42 (-0.64,-0.14) 

-0.63 (-0.78,-0.41) 

0.66 (0.44,0.80) 

0.66 (0.45,0.80) 

Correlation with MSWS (95% Cl) 

-0.50 (-0.70,-0.24) 

-0.42 (-0.64,-0.14) 

0.64 (0.41,0.79) 

0.62 (0.39,0.77) 

t-test p -value 

0.0229 

0.0067 

0.0142 

0.0026 
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